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Abstract 



The matrix model of random surfaces with c — oo has recently been 
solved and found to be identical to a random surface coupled to a q- 
t^J- ■ states Potts model with q = oo. The mean field-like solution exhibits a 

CO , novel type of tree structure. The natural question is, down to which — 

if any — finite values of c and q does this behavior persist? In this work 
we develop, for the Potts model, an expansion in the fluctuations about 
the q — oo mean field solution. In the lowest — cubic — non-trivial order 
Q*^ , in this expansion the corrections to mean field theory can be given 

■ a nice interpretation in terms of structures (trees and "galaxies") of 

, spin clusters. When q drops below a finite q c , the galaxies overwhelm 

O ■ the trees at all temperatures, thus suppressing mean field behavior. 

D \ Thereafter the phase diagram resembles that of the Ising model, q = 2. 

> 

•rH , 

X 1 1 Introduction 



The infinite-state Potts model on random spherical surface can be solved 
exactly Jl|, || by resumming all orders of its low temperature expansion ||] . 
The problem and its solution are of interest for the following reasons. 

• It can be shown to be equivalent to conformally invariant models with 
c = oo (which can be realized, for example, by coupling multiple Ising 
models to the surface). Thus the solution opens a window onto the 
mysterious c > 1 regime, hitherto inaccessible via either Liouville the- 
ory or matrix models. 

• This equivalence to conformal matter holds in spite of a well known 
fact that on a fixed two-dimensional lattice a Potts model with more 
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that four states suffers a first order — rather than a continuous — phase 
transition. Coupling an oo-state Potts model to a random surface 
changes the order of the transition from first to third; in the case of 
q < 4 the order of the transition correspondingly changes from second 
to third. 

• The reason why the transition for q = oo is of third order is that 
it is no longer a spin-ordering transition; rather it is a tree-growing 
transition, where the nodes of the trees are spin clusters (equivalently, 
"baby universes"). The transition in the random surface Potts model 
from spin-ordering (low q) to tree-growing (high q) behavior is a novel 
one, and deserves attention in its own right. 

Some recent work has given support to the mean field theory at q = c = 
oo. Ambj0rn, Durhuus and Jonsson [||, || (see also Q) have shown that if 
a random surface model has susceptibility exponent 7 s t r , its "polymerized" 
version has 7 s t r = 7str/(7str — 1) ; m agreement with the mean field result 
7str = 5 • Their polymerized surfaces of baby universes are suggestive of the 
trees of clusters found in the mean field model]]] Harris and Wheater J?]] 
have studied multiple Ising models, and in the limit c — > oo have identified 
trees as dominant, and found 7 str = \q A recent numerical study (|| of 
multiple g = 2,3, and 4 Potts models shows that spikiness of the surfaces 
increases with increasing central charge, and that 7 s t r may well approach |. 
Finally, a current numerical study of the Potts model for various values of q 
and for surfaces with tadpoles and self-energies || is confirming the scaling 
of the critical temperature T^T 1 = \ \ogq + 0(1) with large q, and, more 
significantly, details of cluster geometry predicted by the mean field theory. 

For convenience, let us introduce the concept of a skeleton graph: given 
a triangulated surface with some spin configuration, its skeleton is a graph 
that represents the spin clusters as simple vertices, and whose edges are the 
links between different spin clusters.^ Each loop of skeleton graph introduces 
a factor 1/q (or 1/c for conformal models). Thus at q = oo (or at c = oo), 
the only skeletons that are present are trees. When we move down to large 

1 Although baby universes (geometrical structures) are not a priori identical to clusters 
(matter structures), the two seem to fuse at large q or c. Note that the trees in |^| are a 
result of the N — > oo limit, while the trees in jj], ^] are a result of q, c — » oo. 

2 Their trees are made of vertices, not clusters, and thus form a subset of trees-of- 
clusters. There is no contradiction: to find clusters at the vertices of the trees, one has to 
take c — > oo and T — » simultaneously, keeping e~ 2 ^ T c fixed. 

3 So be definition, no vertex of a skeleton may be connected to itself. 
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but finite q, the corrections introduce loops; and while the leading terms are 
universal, the corrections in 1/q are not — their details depend on the details 
of the matter model. This is a familiar picture in 1/D expansions. For 
instance in the scalar A</> 4 theory, the upper critical dimension is 4: below 
this, fluctuations or "loops" become essential, and mean field theory breaks 
down. 

It would be very interesting to know if something similar happens for 
random surfaces. One's curiousity is heightened by the fact that Liouville 
theory seems to be silent on the matter. To investigate the large-but-finite 
q region, we develop a systematic expansion of the model (in its equivalent 
one-matrix form) about the q = oo solution. A priori, we would expect 
that the first non-trivial order in this expansion to be quadratic, as in many 
similar cases. In this order we sum all the skeletons that have one or no 
loops, in other words the q° and q' 1 terms in the partition function. But it 
is not hard to show that this q^ 1 correction is insignificant, in that it does 
not modify the critical behavior of the model. 

The next, cubic, order of our approximation is much richer. Here we 
include all skeletons whose clusters have three or fewer external legs that 
lie on loops. We therefore have terms with arbitrarily many loops, meaning 
arbitrarily high powers of q^ 1 . A new structure makes its appearance: the 
clusters are organized into 3 graphs of spherical topology. These "surfaces" 
of clusters will be called galaxies. The trees are still present: any number of 
them may grow out of any cluster. On top of that, a typical skeleton graph 
includes several galaxies, linked also into a tree. 

This complicated system shows complicated critical behavior. For q = 
oo, there is a magnetized low temperature phase where the clusters are 
large, and there are few of them: essentially pure gravity. At some point 
that phase becomes unstable, and the tree of clusters blows up. Above that, 
the clusters are small, but the tree is large. 

At large but finite q, the above picture of critical behavior is unchanged, 
except for the appearance of a new phase at high temperature, that of large 
galaxies. At low temperature we still have the pure gravity phase of large 
clusters; this still goes into an intermediate phase of large trees. But now, at 
an even higher temperature (which is infinite when q is infinite, and which 
decreases with q), the area of the galaxies — the number of clusters they 
contain — diverges. In this phase, the loops become essential: without loops 
there could be no galaxies. The phase of large galaxies is also akin to pure 
gravity. 
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The tree phase becomes narrower and narrower as q decreases, until 
at some q c , it disappears altogether: one goes from the magnetized phase 
directly into the phase of large galaxies. Although mean field theory is 
modified at any finite q, it is completely destroyed below q c . For q < q c 
the phase diagram is altogether reasonable, in that it resembles the known 
phase diagrams of the q < 4 models: there are two pure gravity phases, at 
low and at high temperature, separated by a multicritical point. 

This paper is organized as follows. Section 2 presents the multi- matrix 
version of the Potts model, its low temperature expansion, its effective one- 
matrix versions, and its q = oo solution. In section 3 we introduce an 
expansion about the q = oo solutions and its truncated versions, and discuss 
its geometrical implications. In section 5 we analyze the critical behavior of 
the simplest nontrivial truncation. 



2 Definition of the model 



The g-state Potts model on a spherical random surface will be defined using 
a matrix model (see [|l(| for a review): 

F q (g, C) = - log j jvfa ■ ■ ■ V<p q exp -N TV I ± £ faG^ + 9 £ 4 I 

(1) 

where the <f>i are N x N hermitian matrices and all matrix integrals are to 
be understood as being divided by N 2 and taken in the limit N — > oo. The 
matrix G defines the matter that is to be coupled to the fluctuating surface; 
in the case of the Potts model it is 
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(2) 
(3) 

e- 2 / T . 



The coupling constant ( is related to the matter temperature by £ 
Finally, the normalization /o is defined by requiring F q (0, = 1- 

The partition function F q can be expanded in powers of £. This is the 
so-called low temperature or cluster expansion, as in order Q n we have a 
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sum over all connected configurations of clusters with n inter-cluster links; 
as well as over the internal geometry of each cluster. Defining A = q — 1, 
the first four orders of the expansion are 



F g (g,c) 



1 9 

+ 2 A?r iC + 



+ 7 A7r 2 

4 z 
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Ja 3 tt 2 (^ 2 + 7r n ) 2 + -A 3 7r?(7r 3 + n 12 + 7r m ) + ^A(A - 1)tt| 

Z O 
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(5) 



The picture underneath each term is the corresponding skeleton graph — see 
the Introduction. The various coefficients ir p (p is a partition of an integer 
n, where n is the cluster's external legs) are functions of the cosmological 
constant g that sum over the internal geometries of the clusters. 

To define and to calculate the 7r's, one introduces an external matrix 
integral 

r%,A) = --LTrA 2 -logi-| P0exp-iVTr (\<j> + ^<f> 2 + 5 3 ) (6) 

This integral has been calculated by Kazakov and by Gross and Newman. 
If we expand it in powers of A, we get 



n(<7,A) 



ITi , 1 

vr + ¥ TrA + - 



-^■Tr A + ^(TrA) 



— Tr A 3 + — !§Tr ATr A 2 + — ^ (Tr A) 3 

N N 2 N 3 



+ 



(7) 



The coefficients Tr p (g) are the ones we encounter in the low temperature 
expansion @). In n-th order, we find ir p for all partitions p of n. The 
coefficient ir nin2 ... is a connected Green's function of the one-matrix model 

H(g, 0) with n\ +n 2 H external legs, on each of which we place the matrix 

A; the graphs can be twisted in various ways, which gives rise to the different 
combinations of traces. For example, tt2 sums graphs such as 



while it n sums graphs such as 



(TV A) 5 



The untwisted coefficients 7r n are just the n-point connected BIPZ Green's 
function []TT| , while the twisted coemcients have not been discussed before 
in the literature. They can all be calculated exactly, and the first few are 
given in Appendix A. 

The complicated linear combinations of the coefficients ir p appearing in 
the low temperature expansion @ are due to the large N limit of our model. 
It will be seen below that there is a way to organize these linear combinations 
that will clarify their nature and simplify their calculation. 

Here we briefly summarize the solution of the Potts model in the q — > oo 
limit. As it can easily be seen that the coefficient of £ n in the partition 
function is an n-th degree polynomial in A, the proper scaling for the tem- 
perature is C ~ 1/A. Defining £ = a/A and keeping a fixed0 while letting 
q — > oo, it was found that only skeleton graphs that are trees contribute to 
F QO (g, a). These can be resummed to give 



Foc(g,a) 



lim 



1 



log 



n(s,y) 
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2a 



dx exp —9 



x 
2a 



U{g,x) 



U(g,y)=U(g,yl) 



where the function y(g, a) is a saddle point of F a 
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y 



a d-y U(9 ' y) 



(8) 
(9) 

(10) 



The general features of the q = oo model's critical behavior have already 
been described in the Introduction; see Q for details.^ 

Kazakov [12] has pointed out that the random surface Potts model can 
also be formulated as a one-matrix model with the action 



S n = A 



lTr0 2 -iv(l + i)n((l-a/A) 3 / 2 , 



a/A 
(11) 



4 This is the most interesting way of taking the q —* oo limit. If we keep £ fixed instead, 
for example, we will find only one of the several phases of the model. 

5 This limit can be called a mean field solution because it can be arrived at by (an a 
priori inexact) procedure of substituting (f>i(f>j —+ <&4>i in the multi-matrix model (Jl]) and 
solving self-consistently for additionally, one has to take the (not a priori obvious) 
ansatz that $ is a multiple of unity, as in Wilson's mean field solution of QCD. 
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The simple derivation is given in Appendix B. This effective action will be 
called the "natural" one. For the purposes of the present calculation, it is 
useful to introduce a slightly different, "artificial" effective action 



S a = A 



-pTr d> H Tr 

H Y 2a 



(12) 



The difference between these two effective models, their relative merits and 
drawbacks, are fully explained in Appendix B. 

In going from the multi-matrix model (|l|) to the effective one-matrix 
models (|ll]) and (|l^), we trade the complexity of having several matrices 
for the complexity of the new one-matrix potential — see eq. (0). But what 
is really happening is that the vertices of the effective models are playing a 
completely different role than the vertices of the original model. In the orig- 
inal multi-matrix model, the vertices represented points on the discretized 
worldsheet. In the effective models, the vertices — which are complicated 
functions of the original cosmological constant g — actually represent entire 
spin clusters. This is possible because a twisted cluster is equivalent, from 
the large- ./V fatgraph point of view, to a nonlinear vertex 



such as (TrA) 2 . This is why all possible nonlinear terms appear in the 
expansion (^) of the potential II(<7,A); or, seen from the other direction, 
why clusters of every possible twist appear in the low temperature expansion 

(D-fl 

One could rederive the q, A — > oo limit in a quick and dirty way by 



noticing that when the model (11) is reduced to an eigenvalue problem, the 
charge of the eigenvalues will scale as A -1 , and therefore the width of the 
eigenvalue spectrum will scale as A -1 / 2 . At A = oo, the eigenvalue density 
will be a delta function — precisely the content of equation (||) . 



3 Fluctuations about the mean field 



So at q = oo the eigenvalue density of the effective models (11) and (|12| ) 



collapses to a delta function. When q is large but finite, the eigenvalues 

6 The procedure that led from multi-matrix model to the effective one-matrix model can 
be applied to any — not only Potts — matter. Only the Potts model, however, is symmetric 
with respect to any relabeling of matter states; that is why its effective model is a one- 
matrix model. Effective cluster models for other types of matter are multi-matrix models. 
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spread into a narrowly peaked distribution. The main idea of this work is to 
investigate the shape of that distribution by expanding the effective models 
about the scalar solution 

4> = yl + ^ (13) 

in powers of ip. One type of critical point will result from singularities of 
the equations for the eigenvalue distribution of ip. Once an approximate 
distribution is found in this way, we will plug it back into the full effective 
model to determine critical behavior; another type of critical point will 
result from singularities of the potential II for a given ip. For simplicity, 
the following calculations will be shown only for the artificial model (]i~2|) ; 
the calculations for the natural model (|ll|) differ in trivial ways, and their 
result will be given at the end. 
The action (pi) now reads 



S a = A 



(14) 

The reader is reminded that the q = oo partition function is a sum 
over free trees, and its saddle point y(g, a) is a sum over planted trees: trees 
with a marked vertex of order one. Explicitly, 



y(g,a) = 7ria + 7ri(7r 2 +7rn)a 2 + 



7Tl(7T2 + VTll) 2 + -7T 2 (7T3 + TT\2 + Kill) 



This low temperature expansion of y can be represented graphically as 

°- 1 + 



+ ••• 

(15) 




We are led to consider the expansion of the external matrix integral about 
a scalar: 



n(g,y + i;)=n (g,y) + ^y±Trr 



1 

+ 2 
1 

+ 6 



N 

n2(g ^W+ ni ^W 2 



(16) 



N 

113(9, y) 

N 



N 2 

TrV 3 + ^^TrV-Tr^ 2 + ^^(Tr^f 



N 2 



N 3 



+ 
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The vertex functions n p can be calculated explicitly, and are given in Ap- 
pendix A. To understand what they are counting, we can expand several of 
them in powers of y: 

1 1 

no(0,2/) = 7ro + 7riy+-(7r2 + 7rii)y 2 + -(7r 3 +7ri2 + 7r m )y 3 H 

2 b 

TT I \ I ( . 7I"12\ . 1 / . 7T13 . 7T22 . ^112^ 2 , 

n 2 (ff, y) = vr 2 + I 7T 3 + — I y + - I vr 4 + — + — + — jy +■■■ 
113(9, y) = 7i"3 + (ji"4 + y H , etc. (17) 

The linear combinations of 7r p 's in Ho, 112,113,... are precisely those that 
occur in the low temperature expansion (Q). The ra-th coefficient of Ho 
counts clusters that have n external legs, each of which is connected to a tree, 
i.e., none of its legs reconnect. By contrast, the n-th coefficient of n 2 counts 
clusters that have two legs that lie on a loop (and will eventually reconnect), 
and n other legs that are connected to trees; and so on. Thus a cluster 
with three external legs, for example, may appear in the low temperature 
expansion as ir^, ir^ + 7Ti 2 /3, or tt^ + 7Ti2 + 7Tm, depending on whether it has 
3, 2, or legs that lie on loops. 

The expansion of action (14) in powers of -0 is therefore[] 



% = - N 
A 



" " W^^^ " ^ (TrV;)3 + ' ' ' (18) 



According to the above discussion, the vertex H3 Tr tp 3 is really an infinite 
sum of vertices ( |l~7|) , and can be graphically represented as 

n 3 (y,y)TrV 3 ^ pp + =j= + \ + (19) 

while a non-linear vertex such as II12 Tr ip Tr ?/> 2 may be represented as 



'• 2 ^> 



n 12 ( ff ,y)Tr0Tr^ 

II II 

(20) 



7 Two terms cancel due to (^). 
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The model fllq ) is an expansion about the tree-like, mean field solution 
cf) = y. The geometrical meaning of this statement is reflected in the above 
pictures: trees are automatically included in the vertices of the new matrix 
model. The graphs of this matrix model can be pictured as a dense network 
(in fact, a spherical surface) of clusters, from the vertices of which sprout 
trees. 

To proceed farther with the matrix model ( |l8|) we are obliged to truncate 
the potential H(g, y + ip)- The simplest possibility is to truncate after the 
quadratic term. We take p = in this case, and immediately evaluate the 
partition function^] 

* = PD 

The new 1/A terms in the partition function just count the one-loop con- 
tributions to the low temperature expansion (||). It is clear that they can- 
not effect the critical behavior of the model: a large tree with one loop 
is indistinguishable from just a large tree. This can be demonstrated by 
noticing that the only way in which the 1/A terms can effect the critical 
behavior of the model is through a divergence of the logarithm; but, since 
II2 < c^IIo = II2 + II11 (see Appendix A), and 1 — ac^Ho > in the physical 
region of the q = 00 model (see Q), this is impossible. 

We are therefore obliged to keep higher order terms in the potential 
n(<?, y + ?/>). The meaning of a truncation after ip n terms, for n > 3, is as 
follows. Consider a cluster with A links leading to other clusters. Of those 
links, At are to trees and Xg = A — Xt are to loops; meaning that if one of 
the Xt tree links is cut, the skeleton will be disconnected. For instance, the 
cluster indicated by the arrow 




has Xt = 2 and Xe = 3. So: if we truncate the action (|18|) after the ip n terms, 
we discard any skeleton that has a cluster with Xi > n. Note that Xt does 
not matter: as discussed above, each vertex of (|l8|) automatically includes 
arbitrary numbers of trees. 



We now consider the model (18) in the cubic approximation: we discard 



all terms of order ■0 4 or higher. In this approximation the clusters — which 

8 The last term in (^l]) is there to cancel cluster self-connections; it should have been 
put in by hand as a constant in the action. 
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are themselves spherical </? 3 graphs — are arranged into spherical </? 3 graphs 
that also have trees of clusters growing from their vertices (given by powers 
of y: see (1£), (20)). This structure will be called a "galaxy." A typical 



graph in the cubic model has several galaxies, with the nonlinear terms 
(Tr</>) 2 , Tr Tr </> 2 , and (Tr^>) 3 connecting the galaxies into a tree structure 
(not to be confused with the trees that grow out of the vertices inside the 
galaxies). The cubic approximation differs radically from the quadratic one 
in that it does not simply sum over graphs with one more loop and with one 
higher order in A -1 , but rather it sums over a large class of graphs with an 
arbitrary number of loops. In fact, it sums over all graphs, provided that no 
cluster has more than three legs lying on loops; for example, the following 
are simple skeleton graphs that are not included in the cubic approximation: 




The solution of the cubic model will proceed along the same lines as the 
simpler model solved in [13]. First off, we set 

n3_n 1 2 

y 2a 6a v ' 

to cancel the cluster self-connections. Then we follow the standard pro- 
cedure, reducing the problem to eigenvalues and introducing a continuous 
approximation to the spectrum \{x). Using procedure that is exact in the 
N —* oo limit we introduce the two moments 

c= (jj TripJ =Jdx\(x), d= (^—Tr^J = Jdx\ 2 {x) (23) 

as parameters, for which will we solve self-consistently at the end. The 
saddle point equation for the action reads 

1 I" dy (24) 



6S 

Jx 


= q + m\(x) + 3pA 


V 


n 3 


A 
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m 
~A 


i ^ cn 12 
= — + n 2 + — ^ 

a 6 


Q 
A 


dU 12 

= p + c nn + — H 





w 2 J X(x)-X(y) 

(25) 

(26) 

+ ^ (27) 
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keeping in mind that the lip's are (complicated) functions of g and y{g,a). 
The eigenvalue density for this problem is 



1 , ^ / 4 /, <t \ 2 



„ (A) = _ (m + (T + 3p^— (28) 

where the variable <r satisfies the equation 

18p 2 + a(m + a) (m + 2<r) + 3p<?(m + 2<r) = (29) 



Using the eigenvalue density (|28|) , the self-consistency conditions on c and 
d become 

3P > a (30) 



(m + 2<t) 2 3p 
m + 4ct / <t \ 2 



d = Rw + y (3i) 

The reader is reminded that the eigenvalue density ( p8|) is for the ?/> 
matrix. Since a scales as A , for large A the eigenvalue density for eft in 
the cubic approximation is centered on the q = oo value y with a displace- 
ment cr/3p ~ C(A _1 ). The width of the distribution is y/4/(m + 2a) ~ 
0(A -1 / 2 ), as previously advertised. 



4 Critical behavior 

The cubic approximation to the Potts model is expected to have three critical 
phases: the phase of large clusters, the phase of large trees of clusters, and 
the phase of large galaxies of clusters. 

First, the regime of large clusters. This is the magnetized, low tempera- 
ture phase. As we have an approximate form for the eigenvalue density (ES 



we can plug it back into the full action (12) and look for the singularities of 
the potential H(g, cp). The reader is reminded [Q] that for q = oo, the eigen- 
value density is a delta function at y(g, a); the potential in that case can eas- 
ily be reduced to the "standard" form <j) 2 /2+G<j) 3 , with G = (l-Ugy^^g. 
At the critical point 



G = G C = 1/V 108V3 (32) 

the area of the clusters (but not their number) diverges. The same thing 
happens in the case of the spread out density (|28|), though the computation is 
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somewhat more difficult: we are looking for the critical point of the external 



matrix model (^). The result can be found in [14]. The only dependence on 
the external matrix is through the moments 



d\u(\)(x- \y k/2 (33) 



where the variable x satisfies the implicit equation 

x = ^- - y/3gai(x) (34) 

For the eigenvalue density, these moments can be calculated explicitly (as 
functions of x): 



a k = { x-y-- 



1 k k 

2 + 4' 4'"' {l + 2a)(x-a/3 P y 



2Fi[- + -,-;2; M ; n _ w _ 777^2 ) (35) 



3/cp ^ (\ k k 4 

+ ^ , n wf- 77TT 2-Pl hr + t. 1 + t; 3 ; 



2(1 + 2ff) 2 (x -a/Zp) \2 4' 4' ' (1 + 2a)(x - a/3p) 2 



The criticality condition is the point where the x equation fl34[) ceases to 
have solutions. Algebraically this can be expressed as 

3gal = 4 (36) 

It can easily be checked that when q, A — > oo, and consequently cr/3p, 4/ {m+ 
2a) — > 0, these new criticality conditions (|34]) and (|36]) reduce to (|32|). It 
will be seen in the phase diagrams that as q approaches infinity, the curve of 
large clusters will approach the corresponding q = oo curve, given by solving 
eqs. ( |3^ ) and ( |10| ) simultaneously. 

In the large cluster phase matter is manifestly magnetized, and is there- 
fore decoupled from the geometry; so this is a phase of pure gravity, with 
7str = — \ ■ Another way of saying the same thing: in this phase the number 
of clusters is small, and few clusters is the same as one cluster, as far as 
critical behavior is concerned. 

Second, the phase of large trees of clusters. There are, in principle, two 
ways of getting there. At q = oo, the trees become large when the saddle 
point equation for y(g,a) (|io|) ceases to have a solution §; in other words, 
when 

dy 2 



a^U(g, y)=a [U 2 (g, y) + U u (g, y)} = 1 (37) 
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The other way to get large trees is to make the tree of galaxies (of clusters) 
blow up. This is not exactly the same as a large tree of clusters at q = oo, 
but very nearly the same, as in general the galaxies at the vertices of this 
tree will be small; and as far as critical behavior is concerned, a small galaxy 
is no different from a galaxy of size one, in other words, a single cluster. This 
intuition is supported by the fact that both large tree phases have 7 s t r = +\] 
and as q — > oo, the curves of the two phases approach each other. 

To calculate the location of the second large tree phase, we can look for 
the point in the cubic model (24) where the string susceptibility d 2 F/dp 2 



diverges [O]. Equivalently, and computationally simpler, we can look for 
the divergence of the "pseudo-susceptibility"^ 

dc 

c' S -^oc (38) 



This is most simply done by differentiating the eqs. (29-31) with respect to 
p, and solving the resulting linear system for c', d', and a'. As the result is 
unwieldy, and in any case can easily be derived, it will not be given here. The 
same procedure, if applied to the simple quartic model in Jl3|| , reproduces 
the result for the large tree phase given in that work. 

In the large tree phase one has an ensemble of clusters of different size. 
As one moves toward the magnetized phase, clusters of higher and higher 
area begin to be included. At the multicritical point, some moments of the 
cluster area diverge, a sign of a continuous spin ordering phase transition PI 
The susceptibility exponent in the bulk of the tree phase is 7 str = 5, the 
generic exponent for tree growth. In the tree phase, and only in the tree 
phase, the clusters are probably identical to the geometric structures known 
as "baby universes" ||]. Therefore the transition to large trees from high 
temperature can be interpreted as a magnetization of baby universes (but 
not of the whole surface); and the transition to trees from low temperature 
can be interpreted as a break-up of the magnetized surface into magnetized 
(but in different directions) baby universes Q. 

Finally, the cubic model should also have a phase of large galaxies. This 
occurs simply at the normal, BIPZ, critical point of the model (p4|). This is 

9 d' or a' would do just as well. 

10 For the Ising model on a fixed two-dimensional lattice, for example, the second moment 
of the cluster size distribution is the lowest one to diverge. M] 
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simply given by 



r (--iw/ 4 rc vWi (39) 

This phase is radically different from anything that exists at q = 00; it 
is in this phase that loops become important. Interestingly, it bears some 
resemblance to what must be the case for small q: the spin clusters at high 
temperature do not form trees, but complicated planar networks with many 
loops. 

The large planar networks of clusters in the galaxy phase resemble the 
configurations of a low q model at high temperature. As we go to higher 
temperature, moreover, the clusters get smaller, as matter fluctuations are 
once again decoupled from the geometry. The large galaxy phase is therefore 
pure gravity, once again. The susceptibility exponent is 7 s t r = —5, the 
standard pure gravity value for nonlinear models in this phase. 

We are now ready to study the complicated transcendental equations of 
the cubic approximation. The procedure is as follows. First, fix the three 
adjustable parameters: A = q — 1, g, and a. Then solve the saddlepoint 



equation (10) for y(g,a); this fixes all coefficients in the cubic action (|li 
Next, solve the cubic equation (|30| ) for c as a function of a, plug into ( |29| ) 
and solve numerically for a, taking the root closest to zero. Finally, check 
the criticality conditions: (|3^ ) and (j36|) for the phase of large clusters; (|38|) 
for large trees; or (|39[) for large galaxies. 



The following is a schematic, not-to-scale sketch of the resulting critical 
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curves. 



a 




g 



If drawn to scale, the curves are rather close to each other and hard to 
distinguish by eye. As q — > oo, the curves reassuringly approach the exact 
q = oo critical curve, where the transition from magnetized to tree phases 
occurs at g = 1/V288, a = 1. As we move away from q = oo, the new 
critical regime of large galaxies appears. As q gets lower, the phases of large 
galaxies and large clusters approach each other, eating away at the mean 
field, large tree phase. 

The susceptibility exponent 7 str is — \ in the magnetized and galaxy 
phases, as discussed above, and +\ in the tree phase. At the two multicriti- 
cal points 7 s t r = \- Further, we find a = — 1 at both multicritical points, as 
shown in Appendix C. The equality of the exponents at the two multicritical 
points is not a coincidence: the model at one of the points can be mapped 
onto the model at the other — see the appendix for details. 

At q = q c (q c ~ 133 for the natural model, and q c ~ 108 for the artificial 
model), the phase of large trees ceases to exist. Thereafter, the magnetized 
phase of large clusters goes directly into a phase of large galaxies. The mean 
field behavior is entirely eliminated. Let g c i(q) be the multicritical point 
between phases of large galaxies and large trees, and g C 2(o) the multicritical 
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point between the phases of large trees and large clusters. Ag c = g C 2 — 
g c i is a measure of "how much of mean field theory is left" at a given 
value of q. Here it is plotted for both the natural and the artificial models. 




For q < q c , the two pure gravity phases (magnetized and large galaxies) 
are joined directly onto each other, with a multicritical point in the mid- 
dle. This phase diagram bears a striking resemblance to the known phase 
diagram for, say, q = 2 [16|. 



5 Discussion 



The results obtained for q c cannot be considered as anything more than 
order-of- magnitude estimates. The real result of this calculation is the dis- 
covery of the mechanism for the break-down of mean field theory: the di- 
vergence of galaxies. The resulting critical behavior interpolates smoothly 
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between the somewhat odd phase diagram at q = c = oo and the well-known 
phase diagram for q < 4, c < 1. Another result of this calculation is the like- 
lihood that q c is finite, that mean field behavior holds for some finite Potts 
models (with modifications, as has been shown). 

The two natural questions at this stage are: is the present result merely 
an artifact of the truncation of the action (|l^), or are the higher-order 
terms essential? And does a similar picture hold for conformally symmetric 
models, which, after all, reduce to the same mean field theory at c = oo? 

The first question should not be too hard to answer. It seems that we 
would not find any qualitatively new behavior by including the higher-order 
terms, which would simply add other types of vertices to the galaxies but 
would not modify the essential structure. Mean field theory should always 
be destroyed by the divergence of the galaxies. It is true that quartic or 
higher models will add new multicritical points; but as we only have the 
three parameters q,g, and a to tune, we should not be able to reach them. 
To repeat the present calculation in the quartic approximation should be 
somewhat tedious, but not too difficult. 

The second question is more difficult to answer, and the author does not 
have any suggestions on how to deal with it. The immediate difficulty is that, 
although effective cluster models such as (|TTj) and (|l^) can be introduced 
for any matter model on the surface, they are multi-matrix rather than one- 
matrix models. Still, they closely resemble the c = oo model, so perhaps it 
is the right way to start. 

The author thanks the Niels Bohr Institute for its hospitality, and Jan 
Ambj0rn and Gudmar Thorleifsson for stimulating discussions. Marc Pot- 
ters provided many of the ideas in the early stage of this work. 

Appendix A 

Here we give expressions for the first few Green's functions itp(g) in the 
expansion ([?]) of the external matrix integral II(<7, A), defined in @. Using 
Gross and Newman's notation, the external matrix integral is 



I%A) 



1 




2 



2N 2 



a,b 




(40) 



18 



where the A a are eigenvalues of A, fi a = \Jx — A a , er& = jj J2a( x ~ ^a) 
and the variable x satisfies the implicit equation (|34|). 
First we expand x in powers of A: 

y/5g TrA 

95 (s/Zg - ^T) (TrA) 2 + 3^ TrA 2 + 



4x (V¥-2xf) 3 N2 4, (V%-2xf) N 
where xq satisfies the cubic equation 



x 



1 

125 " 

1 1 / r- , 26 26 

1 V 3 sin cos — 

18s 36 V 3 3 



tan = W — -1, go = , = 

WW ^108^ 

= T ^-6g-216g 3 ---- 
12p 

Plugging the series for x back into (flOl), after some algebra we find 



^ fxo _L l_ 

" ' \l 3g 4x 6g 



3g + 108# 3 + 7776/ + 



7T 2 



7TH 



7T 3 



2 



2^/3gxo 16x, 
27/ + 1944/ + 



16xq 



9/ + 1296/ + • • • 



65 + 540/ + 58320/ + • • • 



16xq 
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324g 3 + 69984/ + • 



7T111 



3<7 



8x3(4x[] - 3g) 
7776/ + • • • 



+ 



35 



4xf (4xg 



3 5 ) 



The corresponding coefficients H p in the expansion of H(g, y + A) (|16| ) 
can be obtained in a similar manner. The above expressions for ir p (g) go 
directly into H p (g, y) with the substitution xo — > xo — y, while xo itself now 
satisfies 

•n, ±- (44) 



^0 - y 



and is given by 



x 



1 y 1 

h - + 

18s 3 



tan 9 



1, G 




26 



cos ■ 



(1 - 12 5 y)3/4 



(45) 



Appendix B 

Consider a model with an extra matrix d> and the action 



+ 



(46) 



with a ? , 6g given in (||). On one hand, the extra matrix can be integrated out 
to give exactly the model ([[]); one the other hand, the fa can be integrated 
out in (46) to give a one-matrix model with the action 



Sn = A 



— Tr 6 2 - N 
2a 



i + ^)n((i 



a/Af/'g 



1 - a/A 

(47) 

Whereas the vertices of the original multi-matrix model ([l]) are the points of 
a discretized worldsheet, the vertices of the effective one-matrix model (47) 
are spin clusters. 

The partition function of the effective one-matrix model (|47|) can be ex- 
panded in powers of a (not a) to reproduce the low temperature expansion 
given above (||). The reason that this works is that a vertex like (Tr</>) 2 , 
such as occurs in the action fl47|), is equivalent — as far as large N behavior 
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is concerned — to a twisted cluster, as explained in the text. The correspon- 
dence between the two expansions is still not completely trivial, because in 
the low temperature expansion (^|) 

• the propagator between clusters is a, rather than d; 

• a cluster is never linked to itself — by definition; 

• the coefficient of a skeleton graph is a polynomial in A" 1 , whose lead- 
ing term is A*' - ^ -1 (where v and £ are the numbers of vertices and links 
in the graph)0; if there are no loops of length greater than two present, 
the polynomial is just a monomial, but otherwise there are subleading 
terms that arise due to the finite volume of target space — for instance, 

the coefficient of is A _1 (l - A -1 ). 

Evidently, the 1/A terms in the action ( f47| ) — the correction to the coefficient 
of II and the renormalization g — ► (1 — a/A) 3 / 2 g — conspire to compensate 
for the above three discrepancies. 

In the calculation to be performed, the potential II in the effective one- 
matrix model (|47|) will be truncated after the cubic term (only the terms 
written out in equation (0) are included). It is instructive to consider the 
effect of this truncation on the low temperature expansion of the effective 
model. Of course the clusters will now have no more than three neighbors. 
Furthermore, consider, say, the coefficient of a 2 : 

r^+^+^-i^^+f) (48) 

The first two terms reproduce the low temperature expansion (Q). The last 
term, however, compensates for the tadpole 

CXD 

that would have been generated by the t/> 4 terms. This over-compensation for 
tadpoles is somewhat undesirable. We can consider an alternate, "artificial" 
action (as opposed to the "natural" action (|47|)) 



S a = A 



-p Tr cf> + — Tr 2 - NU(g, 
2a 



(49) 



11 Assuming that the expansion (M is rewritten in terms of a, rather than c 
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with p adjusted to cancel the simple tadpoles (i.e., the cluster self-connec- 
tions), in this case to — (iri + ^3/2 + 71-12/6)0. Strictly speaking, a constant 
term should also be added to S a to cancel all the cluster self-connections; but 
as this term will have no effect whatsoever on the calculations here, it will 
be omitted. The partition function of the artificial model, when expanded 
in low temperature series, reproduces the low temperature expansion (||) 
exactly, except that the coefficient of each graph is precisely A v ~ e ~ 1 , so that 

the coefficient of / \^ > f° r example, is A -1 rather than the correct A~ 1 (l — 
A -1 ). The artificial model (49) does not know about the target space of the 
original Potts model, and so it cannot reproduce its excluded volume effects. 
On the positive side, the potential II in the artificial action can be truncated 
at any point without introducing tadpole overcompensation terms such as 

To summarize, there are two effective one-matrix models for the Potts 
model, the "natural" model (|47|) and the closely related "artificial" model 
(|49|). If the potential H(4>) in these models is truncated, slight discrepan- 
cies (different for the two models) are introduced into their low temperature 
expansions, as compared to the series These discrepancies are not im- 
portant when q is large. For instance, the critical behavior of the natural 
and artificial models is identical, up to slightly different critical parameter 
values. 



Appendix C 

Here we sketch the calculation of the critical exponent a. It is defined by 
the singular behavior of the partition function F ~ \g — g c \ 2 ~ a , where we 
constrain the approach to the multicritical points to lie on a critical curve. 

At q = 00, we can calculate a exactly for the multicritical point between 
the large tree phase and the magnetized pure gravity phase. The partition 
function is given by eqs. ([|) and (10). The multicritical point occurs at 



g* = 1/v 288, a* = 1. We vary g, and constrain a to lie on the critical curve. 
When g > g*, this is given by 

a+i9) ~ l-V3(9/ 9 o)y* + ( 9 /9o)^ m 
where g$ = l/108\/3 is the pure gravity critical point; when g < g*, the 
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critical curve is given by the equation ((37 



d 2 

a_—n(g,y)=a-[n 2 (g,y)+n u (g,y)] = l (51) 



in conjunction with eq. (P). The partition function turns out to be analytic 
on either side of the multicritical point. Putting Ag = g — g*, and using the 
expressions given in Appendix A, we find 

F + {g) = hm F^g, a+(g)) = {^~\ lo % 2 ) + 96A ^ " 2432V2A 5 3 + • • • 

F..(g) = lim F^g^a^g)) = (^r-\ log 2] + 96A<? 2 - 5760^2A 5 3 + • • • 
9/9* \24 2 J 

The discontinuity in the third derivatives means that the transition is of third 
order, with a = — 1. The corresponding transition from the magnetized 
to the tree phase for 00 < q < q c is similar, and should have the same 
exponents. 

At the other multicritical point, the boundary between the tree and 
galaxy phases, the geometric fluctuations of the clusters (as given by the 
coefficients H p in action ([HI)) no longer matter; what drives the transition 
is the inter-cluster dynamics of the galaxies, as well as the inter-galaxy dy- 
namics. To calculate the critical exponent we therefore may ignore the com- 
plexities of the coefficients H p (g,y) in the action, and consider any simple 
non-linear matrix model which displays similar behavior. As a first example, 
take the quartic model 



S = TV (^ + ^j + ^(Tr</>T (52) 

that was solved in [O]; the exponent a has never been given in the lit- 
erature.^ The phase of "large galaxies" (7 str = —\) occurs on the curve 
g' = </_(<?) and the phase of large trees (7 str = +i) occurs on the curve 
9' = 9+(g), with 



fJ ~ ' 32 

3 



g - 1 (l + v/T+64^) (53) 



9- = -9<?-7V=% (54) 



12 In the context of the original work frL3|] , this exponent is meaningless, as there is no 
outright matter coupled to the surface. 
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The multicritical point occurs at = — 3/256, g'^ = —9/256. As above, we 
fix g' = g'± and find on either side of the multicritical point 

/29 1 , 8\ 64 . 8192 . 2 8388608 . , 

F+ = log Ao H Aq 2 H Af H 

+ V72 2 s 3y 27 y 27 y 243 y 

/29 1 , 8\ 64 A 8192 2 8388608 . - 

F_ = log Ag H Ao 2 An 3 H 

V72 2 6 3/ 27 y 27 y 243 y 

Therefore the transition is again of third order, with a = — 1. 
As another example, consider the nonlinear cubic model 

5 = Tr Q^ 2 +# 3 ) + |^(Tr<A) 2 (55) 

which is closer to our effective model, eq. (|i"8[). Of course it has exactly the 
same critical behavior as the quartic model (^), which can be shown in an 
identical manner. It is more instructive, however, to consider the partition 
function of the model ( |55| ) directly: it is simply given by 

F = U(g,cg')-^- (56) 

where c is the one-point function, given by the self-consistency condition 

c=/^Tr^)=U 1 (g,cg') (57) 

exactly analogous to (|30|). If we put y = eg', we obtain precisely the partition 
function of the q = 00 model @, along with its saddle point condition, eq. 
(|j~0|). Therefore a = — 1, as calculated above. 

We have therefore mapped one of the multicritical points onto the other. 
A similar mapping holds for the full action (p!8|), giving a = — 1 at the 
transition between the tree and galaxy phases. This mapping points up the 
similarity of the magnetized phase to the phase of large galaxies: in the 
former the individual clusters blow up, while in the latter it is galaxies of 
clusters that diverge. Both phases are, essentially, pure gravity. 
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